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Why does the environmental influence on group and cluster galaxies 
extend beyond the virial radius? 
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ABSTRACT 

In the local Universe, galaxies in groups and clusters contain less gas and are less likely to be 
forming stars than their field counterparts. This effect is not limited to the central group/cluster 
regions, but is shown by recent observations to persist out to several virial radii. To gain 
insight into the extent and cause of this large-scale environmental influence, we use a suite of 
high-resolution cosmological hydrodynamic simulations to analyse galaxies around simulated 
groups and clusters of a wide range of mass (logio Mhost/Mf, = [13.0, 15.2]). In qualitative 
agreement with the observations, we find a systematic depletion of both hot and cold gas 
and a decline in the star forming fraction of galaxies as far out as ^ 5 r2oo from the host 
centre. While a substantial fraction of these galaxies are on highly elliptical orbits and are 
not infalling for the first time (~ 50 per cent at 2 r2oo, independent of host mass) or are 
affected by 'pre-processing' (~ 20 per cent of galaxies around groups, increasing to ~ 50 
per cent around a massive cluster), even a combination of these indirect mechanisms does not 
fully account for the environmental influence, particularly in the case of the hot gas content. 
Direct interaction with an extended gas 'halo' surrounding groups and clusters is shown to be 
sufficiently strong to strip the hot gas atmospheres of infalling galaxies out to ^ 5 r2oo- We 
show that this influence is highly anisotropic, with ram pressure along filaments enhanced by 
up to a factor of 100 despite significant co-flow of gas and galaxies. 

Key words: galaxies: clusters: general — galaxies: evolution — galaxies: haloes — galaxies: 
interactions — galaxies: intergalactic medium — galaxies: ISM 
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1 INTRODUCTION 

There is strong observational evidence that the internal properties 
of galaxies depend on their local environmen t. One well-kn own 
example is the morphology-density relation ( lDresslei| [T980l) by 
which early-type galaxies are more common in high-density en- 
vironments such as the central regions of clusters, whereas late- 
type galaxies dominate the field population. While much of this 
apparent trend is explained by galaxies in groups and clusters be- 
ing typically more massive than field galaxies, combined with a 
correlation between galaxy mass and morphology in the field, there 
is still a significant difference in the morphologies of galaxies of 
fixed stell ar mass in the field and those in massive grou ps and clus- 
ters (e.g.. iKauffmann et alJllOOi : iBlanton et alj|2005h . A similar 
relation holds for galaxy colours: those in denser regio ns are pref- 
erenti ally redder than their isolated counterparts (e.g.. lHogg et alj 
l2004h . As the colour of a galaxy is primarily influenced by its re- 
cent star formation activity, this 'colour-density relation' indicates 
a reduced level of star formation in group and cluster galaxies com- 
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pared to those of similar ste llar mass in the field (e.g.. lBalogh et alj 
l2004|Poggianti et al.ll2006l) . 

In principle, both of these trends may either be due to in- 
trinsic differences between field and group/cluster galaxies (i.e., 
galaxies form differently in close proximity to a massive cluster, 
or a large-scale overdensity destined to become a cluster) or they 
may be the result of a transformation of late-type, star-forming, 
blue field galaxies into early-type, passive, red ones after they are 
accreted by a group or cluster. There are several mechanims that 
could trigger such transformations. The tidal field of the group or 
cluster, as well as interactions with other galaxies, which are nat- 
urally more common in dense enviro nments may strip, r e-shape, 
or even totally disrupt a galaxy (e.g.. iMoore et al.lll999l) . At the 
same time, the high velocity of a galaxy relative to the intra- 
groupZ-cluster medium (icmCI gives rise to ram pressure which 
can remove its cold gas reservoir (so-called ram pre ssure stripping, 
iGunn & Got3ll972l ; lAbadi. Moore & Boweilll999l) and therefore 

' We use the term 'ICM' for hot gas in and around groups/clusters which 
is not bound to infalling galaxies or sub-groups. This includes gas beyond 
the virial radius. 
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shut down star formation. The hot gas haloes surrounding infalling 
galaxies are even easier to remove since these are much more ex- 
tend ed, less dense and th erefore less tightly b ound to the galaxy 
(e.g.jLarson. Tinsley & Caldwe ll 1980: Balo gh. Navarro & Morris! 



I2OOOI : [McCarthy et al.ll2008l) . While not impacting star formation 



directly, the removal of hot gas does put an end to the replenish- 
ment (through cooling) of the cold gas used up in it. The result is 
a delayed decrease in star formation as the remaining cold gas is 
consumed ('strangulation' or 'starvation'). The common theme of 
these mechanisms is that the environmental influence on galaxies 
falling into a group or cluster decreases with increasing distance 
from the centre, as the density of both the ICM and galaxies, as 
well as their orbital velocities, decrease. It is less clear, however, at 
which point during a galaxy's infall one should expect these mech- 
anisms to first 'switch on'. This depends on the detailed structural 
properties of the infalling galaxies as well as that of the host groups 
and clusters into which they are falling. 

There is ample evidence from both observations and simula- 
tions that groups and clusters have no sharp 'edge' to mark the tran- 
sition to the field environment. Instead, their outer regions blend 
smoothly into the surrounding large-scale structure. There is there- 
fore no obvious starting point for the above-mentioned mecha- 
nisms to begin acting on galaxies. A commonly employed bound- 
ary radius for a halo is the 'virial radius' (hereafter ^200)^ which 
is often computed as the radius inside which the average density 
equals 200 times the critical density. Roughly speaking, this ra- 
dius corresponds to the extent of the virialised region of a halo 
in cosmological simulations. There is, however, mounting obser- 



distances significantly greater than r^on (e.g.. iBalogh et alJ|l99S : 


Hansen et alJl2009VHaines et al." 2009': von der Linden et alJl20IC ; 


Luetal, 2012 


; .Wetzel, Tinker & Conroy ,2012: .Rasmussen et al] 


2OI2I: see also 


Bahe et alj|2012h. If environmentally-induced trans- 



formations are indeed responsible for the observed trends, then 
galaxies are evidently affected well beyond the virial radius. 

One commonly identified indirect way in which galaxies can 
be environmentally affected at large distances from the centre of a 
galaxy cluster is through 'pre-processin g' in infalling groups (e.g., 
iBerrier et ai] |2009l : lMcGee et aljlloogh . A second indirect mecha- 
nism is what we refer to as 'overshooting': a non-negligi ble frac- 
tion o f infalling galaxies are on highly elliptical orbits (e.g.. lBensonl 
l2005h that bring them well within r2oo on fi rst pericentric passage 
but back out beyond this ra dius later on (e.g.. lGill, Knebe & GibsonI 
l2005l : lLudlow et al.ll2009l) . The existence of environmental trends 
out to radii well beyond 1200 is therefore not necessarily incompat- 
ible with direct environmental influence being confined to smaller 
scales. 

Of course, the existence of these indirect mechanisms does not 
rule out the possibility that there is also direct environmental influ- 
ence of the group or cluster at distances beyond the virial radius. 
As groups and clusters blend into the large-scale surrounding envi- 
ronment, a galaxy is surro unded by gas even at large distances from 
the host centre (see, e.g., iFrenk et alj [l999l who show that the hot 
gas haloes of simulated massive galaxy clusters extend out as far as 
~ 10 Mpc) This gas will exert a ram pressure force on any galaxy 
moving relative to it, which may be sufficient to remove (some of) 
its gas. 

It is therefore conceivable that all three of these mechanisms 
are operating in concert to reduce the star formation activity of 
galaxies. However, separating them from each other would be a 
formidable challenge for current observations, as there is in gen- 
eral only limited information about a galaxy's velocity and that of 



its surrounding gas. Furthermore, it is extremely difficult to eve n 
detect the hot gas haloes of the nearest galaxies teregmanllioOTh . 
let alone measure their structural properties around more distant 
galaxies falling into groups and clusters. 

Numerical simulations are a potentially promising tool to gain 
further insight, as both velocities and galaxy orbital histories are 
readily available. The demands on these simulations, however, are 
considerable. Cosmological initial conditions are required for real- 
istic galaxies and hosts, and to give meaningful information on ram 
pressure, baryons have to be included in the simulations directly, 
along with realistic physical prescriptions for relevant processes 
that affect the baryons, including radiative cooling, star formation, 
chemodynamics, and supernova feedback. The resolution must be 
high enough to resolve individual galaxies, while the simulation 
must also include rare objects such as massive galaxy clusters. 
High-resolution cosmological hydrodynamic simulations of large 
volumes would satisfy these requirements but their computational 
cost is currently prohibitively high. A promising compromise is the 
use of simulations with 'zoomed' initial conditions, where a small 
region of a large, lo w-resolution simulation box is re- simulated at 
high resolution (e.g.. lTormen, Bouchet &Whitelll997h . 

In this paper, we use a set of simulations following this philos- 
ophy, the Galaxies-In tergalactic M edium Interaction 
Caculations (GIMIC: ICrain et al.ll2009h . By re-simulating se- 
lected regions of the Millennium Simulation jSpringel et aTll2005l: 
which includes dark matter only), chosen to encompass a wide 
range of large-scale environments they include both sparse voids 
and a massive galaxy cluster, as well as numerous less massive 
clusters and galaxy groups. A particular advantage of GIMIC is that, 
over a relatively wide range of stellar mass, the simulated /leW disc 
galaxies have already been shown to have properties in good agree- 
ment with a variety of observational data, including the relations 
between stellar mass and rotation velocity, size, and star forma- 
tion effici£ncy_(defined_as the ratio of stellar mass to total mass: 
see [McCarthy et"d]|20I2d) . The simulations also reproduce the 
observed scalings of hot gas X-ray luminosity with K-band lumi - 
nosity, star-formation rate and rotation velocity jCrain et al]|2010t) . 
as well as the properties of stellar haloes around Milk y Way-mass 
disc galaxies jpont et al.f201 iMMcCarthv et alj20I2al) . With a rea- 
sonably realistic population of field galaxies, these simulations are 
therefore suitable to investigate the processes acting on them upon 
infall into a group or cluster 

This paper is structured as follows. In Section 2, we briefly 
describe the simulations and our method for identifying and trac- 
ing galaxies within them. The extent of environmental influence 
on our simulated galaxies beyond r2oo is shown in Section 3, fol- 
lowed by an in-depth analysis of the underlying physical mech- 
anisms in Section 4. In Section 5 we investigate the influence 
of filaments, before presenting our conclusions in Section 6. All 
masses and distances are given in physical units unless other- 
wise specified. A flat ACDM cosmology with Hubble parameter 



h = Ho/( 100 km s" 'Mpc- 



; 0.73, dark energy density parame- 



ter Slj^ = 0.75 (dark energy equation of state parameter w = — 1), 
and matter density parameter Slyi = 0.25 is used throughout this 
paper. 



2 SIMULATIONS AND ANALYSIS 
2.1 GIMIC simulations 

This work is based on the Galaxies-Intergalactic MEDIUM 
Interaction Calculation suite of simulations (gimic). The 
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reader is referred to lCrain et alj l l2009l see also lSchave et alfcoioh 
for a full description of these simulations; here we only summarise 
their main features relevant to this study. 

The GIMIC simulations are a set of five re-simulations of 
nearly spherical regions of varying mean densi ty extracted from 
the Millennium Simulation jSpringel et alj|200^ . The regions are 
chosen so that at z = 1.5 their average densities differ from the 
cosmic mean by (-2, -1, 0, +1, +2) a, where a is the rms mass 
fluctuation on a scale of 18/i^' Mpc at this redshift. In this way, 
GIMIC includes rare objects such as a sparse void and, of particular 
importance here, many groups and clusters of galaxies, including a 
particularly massive one, with logjg (M200 / Mf^) ~ 15.2 at z = 
at the center of the +2(7 sphere. 

The simulations were carried out at 3 different resolutions: 
'low', 'intermediate', and 'high'. The 'low' resolution is the same 
as in the original Millennium Simulation while the 'intermedi- 
ate' and 'high' resolution simulations have 8 and 64 times better 
mass resolution, respectively. As only the —2a and 0(T regions 
have been run at high resolution (owing to prohibitive compu- 
tational expense), we use the intermediate-resolution simulations 
here. These simulations have a baryon particle mass resolution 



of 

yyissis ^ 1.16 X 10^ /i ^Ai:^, with a gravitational softening that is 
kpc in physical space at z < 3 and is fixed in comoving 

ax- 



1 hr 

space at higher redshifts. Thus, even relatively low-mass 
ies (M* ~ a few lO' Mq) are resolved into several hundred parti- 
cles, making GIMIC suitable to study the interaction between galax- 
ies and groups/ c lusters of a wide range of masses. We note that 
[McCarthy etai] ( l2012bl) have shown that the star formation effi- 
ciencies and sizes of the simulated galaxies in GIMIC are approx- 
imately converged when there are (at least) several hundred star 
particles present (although a larger number is required before the 
z = specific star for mation rates converge, s ee Section 3 for fur- 
ther discussion), while lMcCarthv et alj l l2008h have shown that the 
stripping of hot gas is converged when there are a similar number 
of hot gas particles present initially. 

The simulations were carried ou t with the Tre ePM-SPH 
code GADGET- 3 (last described in ISpringell l2005h and in- 
clude significantly modified prescriptions for star formation 
jSchave & Dall a 'Vecchia'2008'). metal-dependent radiative cooling 
in the presence o f a Haardt & Madau ( 2001) UV/X-ray background 
jWiersma. Schave & Smithll2009l) , feedback and mass transport b y 
Type la and Type II supernovae jPalla Vecchia & Schave|]2008h. 



as we ll as stellar evolution and chemodynamics 1 Wiersma et alj 
I2OO9I) . However, they do not include a prescription for feedback 
due to accreting supermassive black holes, so that massive galax- 
ies withlog|fi(Mt /Ma ) ^, 10.7 in GIMIC suffer fro m 'over-cooling' 
(see lCrainetai]|2009l and lMcCarthv et al.ll2012bl) . As a result, we 
do not expect realistic predictions for massive galaxies and there- 
fore limit our analysis to the range logi()(M*/M0) = [9.0, 11.0]. 
We note that the neglect of AGN feedback in these simulations also 
means that the central regions of the simulated groups and clusters 
will not be realistic in terms of the ICM properties or the size/mass 
of the central brightest galaxy jMcCarthv et alj|201(ll) . However, at 
very large distances from the group/cluster centre, the region on 
which we focus in this study, the cooling time of the ICM is much 
longer than a Hubble time, so that this does not affect the validity 
of our results. 



2.2 Selection and tracing of iiost groups and clusters 

Host groups and clusters of galaxies were identified at redshift 
z = using a standard Friends-of-Friends (FoF) algorithm with 
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Figure 1. Comparison between mass within r2oo ('M200') arid total gravita- 
tionally self-bound mass ('Mbo„„j') of the host haloes used in this work. For 
most systems, there is very close agreement between these two, so that the 
host would be included in our sample with either definition (blue points). In 
general, Mtound is slightly lai'ger than M200 due to bound structures extend- 
ing beyond r2oo, so a small number of systems is only included when using 
the former criterion (green points). 



a linking length of b = 0.2 times the mean inter-particle separa- 
tion. We select as hosts all FoF groups with Mbounj > lO'"^" Mq, 
where Mbound is the mass of all gravitationally bound particles 
within this FoF gr oup as identified by the SUBFIND algorithm of 
iDolag etalJ(l2009h. T his ver sion extends the standard implementa- 
tion of lspringel et al.l ( I2OOII) by including baryonic particles in the 
identification of bound substructures and also allows one to distin- 
guish substructures which are located within still larger substruc- 
tures (i.e., sub-subhaloes, sub-sub-subhaloes etc.) from those which 
are associated with the main subhalo of a FoF group. 

This type of host mass threshold is somewhat different from 
the more commonly used spherical overdensity mass M2oo- The 
reason for our choice is that we have identified many instances 
where the FoF algorithm will link together, e.g., a cluster-mass 
halo and a nearby infalling group-mass halo (see Section |4j2). In 
this case, no value of M200 would be computed for the group-mass 
halo, since it is part of the overall FoF group, while Mbound is still 
well-defined. 

In Fig. [T] we present a comparison between Mbound and M200 
for our host FoF haloes. Most of them trace a narrow sequence with 
Mboiind ~ M200 and so are selected irrespective of which type of 
mass cut is applied (blue points). The former is slightly higher on 
average, due to bound structures extending beyond r2oo> so there 
is a small set of galaxy-group scale haloes that are only included 
in our sample adopting the Mbounj cut (green points). In total, 
our sample includes '--^ 100 host systems with masses in the range 
13.0 <logioMbou„d/ Me < 15.2. 

To identify the progenitors of the host haloes in previous snap- 
shots, we use all the gravitationally bound dark matter particles in 
each FoF group in our z = host sample as a tracer population. 
Using their unique particle IDs we identify the FoF halo to which 
the majority of the tracer particles belonged in the previous snap- 
shot, which we designate as the progenitor. We repeat this tracing 
procedure back to redshift z = 10. 
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2.3 Galaxy identification and selection 

Having identified the host haloes and their progenitors in each sim- 
ulation, we next select 'galaxies' (i.e., self-bound subhaloes) in a 
similar way. Starting at the highest redshift (z = 10), we identify 
for each galaxy detected by SUBFIND its constituent dark matter 
and star particles. In the subsequent ('target') snapshot, we then 
search for the galaxy containing most of the mass of these particles 
and identify it as the original galaxy's descendent. Any subhalo in 
the target snapshot which is not identified as a descendent is taken 
as the starting point for a new galaxy and the process is repeated 
until reaching the snapshot at z = 0. 

In the case that two or more galaxies have the same descendent 
in the target snapshot, we continue tracing only the one contributing 
the most mass; all others are marked as 'accreted' onto this galaxy 
and are not traced further so that no galaxy is counted twice. 

We also take into account the possibility that a galaxy may 
temporarily not be identified as gravitationally self -bound by SUB- 
FIND (e.g.. [Muldrew, Pearce & Power 201 Ij. If a galaxy in snap- 
shot i has no descendent in snapshot ; + 1 , or an identified descen- 
dent accounts for less than 50 per cent of its mass, we repeat our 
seach in snapshot ; + 2 and, provided a descendent is identified in 
this snapshot, we continue tracing the galaxy from there. This is 
particularly important for galaxies moving through the central re- 
gions of galaxy clusters where the high background density makes 
an erroneous non-detection more likely. 

To ensure that we can trace a galaxy falling into a group or 
cluster for as long as possible, we exclude all particles belonging to 
a host when identifying the descendent galaxies. In this way, even if 
the vast majority of the galaxy's particles has been stripped, we will 
then still identify the subhalo made up from the remaining bound 
particles as its descendent, and not the host halo, which would only 
be the appropriate choice if the galaxy had been totally disrupted (in 
practical terms this means that there are less than 20 bound particles 
remaining, at which point SUBFIND no longer classifies a structure 
as a self-bound subhalo). 

Finally, we select for analysis those galaxies whose total 
(bound) stellar mass in at least one snapshot falls within the range 
logio(M, / M0) = [9.0, 11.0]. This results in a final sample of 
~ 30000 unique galaxies over all redshifts. From these, we create 
two sub-sets: our main sample of galaxies in the vicinity of hosts 
(the 'infair sample) is formed by those that are found within 5 r2oo 
from the centre of a host in at least one snapshot; there are 15 000 
galaxies in this set. For comparison purposes, we also form a sam- 
ple of 'field' galaxies, defined as centrals which never come closer 
than 5 1200 to ™ FoF group with Mbound ^ lO'"' M,;,. 

Our reason for identifying hosts at z = and then tracing them 
backwards in time, while tracing galaxies forwards from z = 10 
without any prior selection, is two-fold: on the one hand, being in- 
terested in the influence on galaxies by their host environment, we 
want to select only those hosts that are themselves evolving rela- 
tively undisturbed, without being accreted onto other, more mas- 
sive hosts. Identification at z = and then tracing backwards these 
'surviving hosts' satisfies this purpose. For galaxies, on the other 
hand, we also want to include those identified at z > which have 
been disrupted or merged later on, and take into account the fact 
that the stellar mass of a galaxy may vary significantly over cosmic 
timescales. For this reason, we have chosen to trace them forward 
in time, and only select the actual 'galaxy' subhaloes afterwards. 



3 LARGE-SCALE ENVIRONMENTAL TRENDS AT Z ~ 

Within our traced sample of galaxies and hosts as described above, 
we now look at the extent of environmental influence on satellite 
galaxies. For the purposes of the present study we focus on the 
fraction of galaxies which are star forming and have (hot or cold) 
gas mass fractions exceeding some threshold value. We intend to 
explore a wider range of properties, including colours and various 
measures of morphology, in a future study. 

In Fig.|2]we show the fraction of galaxies in which the ratio of 
total mass of gravitationally bound hot and cold gas to stellar mass 
(i.e., the hot and cold gas mass fraction, respectively) lies above 
a threshold value of 0.1, as well as the fraction of galaxies with a 
specific star formation rate (sSFR) above 10^ " yr^ ' . Both of these 
thresholds are set by the resolution of our simulations and the desire 
to have a fixed (specific) threshold across the entire range of stel- 
lar masses that we explore. Conveniently, the sSFR threshold that 
we adopt is very similar to that which is employed in many obser- 
vational studies: these show a well-defined star-forming sequence, 
which is i solated from passive galaxies b y an sSFR cut at 10^" 
yr"' (e.g.. I Wetzel. Tinker & Conrovll2012l) . To distinguish hot and 
cold gas, we adopt a threshold temperature of T = 2 x 10^ K; for 
cold gas we additionally require a density n > 0.01 cm^-^. Such a 
cut in density and temperature for the cold gas roughly mimics a 
selection of atomic (HI) gas. 

For comparison, we also show the corresponding fractions in 
our field sample as yellow dashed lines. To obtain statistically ro- 
bust results despite the multiple splits imposed on our galaxy sam- 
ple (M^, Mjiost. 1/1200)5 we make use of the fact that each galaxy 
in our simulations is 'observed' in more than one snapshot, but at 
different points during its infall into the host group or cluster. We 
therefore stack the results for the redshift range < z < 0.5 to give 
a total of ~ 50000 data points. We have explicitly verified that us- 
ing a smaller redshift range has no significant effect other than to 
increase the statistical uncertainties. 

The first clear influence shown in Fig.[2]is that of galaxy stel- 
lar mass. More massive galaxies (in the right panels) are overall 
considerably more likely to contain hot and cold gas and to convert 
this gas into stars. In the case of hot gas, this comes as no surprise, 
because more massive galaxies with their deeper potential wells 
can be expected to accumulat e and shock-heat more gas than their 
lower-mass counterparts (e.g.. lWhite & Fren3ll99ll) . That a larger 
fraction of massive galaxies have substantial cold gas fractions and 
are star forming compared to the lowest (stellar) mass galaxies 
in our sample, however, appears to be at odds with observations, 
which show that it is in fact low- mass galaxies that t end to have 
large r cold gas mass fra ctions (e.g., lDutton et al.ll201ll) and sSFRs 
(e.g.. ISalim et al.ir2007h . At lo w masses, this di s crepanc y is likely 
a numerical effect: as shown in lMcCarthv et alj l l2012bh . the z = 
sSFRs in GIMIC are only numerically well-converged for galaxies 
with logio(M* / Mq) > 9.7. This is likely due to the fact that a small 
number of particles results in a poor sampling of the gas density 
PDF, keeping in mind that only the highest-density gas forms stars 
and is classified as 'cold' here. On the othe r hand, the absence of a 
'mass-quenching' effect jPeng et ai]|2010l) in more massive galax- 
ies in our sample is a result of incomplete sub-grid physics (such 
as the lack of AGN feedback). While this is clearly not desirable 
and makes it difficult to draw quantitative conclusions, we can still 
analyse trends with environment to gain qualitative insight into the 
relevant physical processes. 

Leaving aside the dependence on galaxy mass by comparing 
galaxies with similar M* (i.e., comparing trends in the same col- 
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umn in Fig.[2]l, there is also clearly a strong influence of galaxy 
environment. Within each panel, the fraction of galaxies above the 
threshold gas mass or star formation rate increases systematically 
with increasing galaxy distance from the host centre. The behaviour 
is similar for all three quantities under consideration (hot gas, cold 
gas and star formation), with hot gas being affected more strongly 
and out to larger cluster-centric radii. 

A second influence is that of host mass: almost universally, 
galaxies around a massive cluster (red lines) are more depleted than 
galaxies near a low-mass group (black) at the same host-centric dis- 
tance in units of r2oo> although the difference between different host 
masses at the same radius is generally smaller than the radial vari- 
ation. There is also a correlation between the influence of galaxy 
and host mass: While low-mass galaxies (left-most column) exhibit 
strong radial variations in both their hot and cold gas content across 
the whole range of host masses under consideration here, massive 
galaxies (right-most column) show only a very mild effect on their 
cold gas content and ability to form stars in environments other 
than a massive cluster. All of these trends are broadly consistent 
with what is expected for ram pressure and tidal stripping. 

Strictly speaking, the existence of these trends alone does not 
guarantee that galaxies are actually changing during infall: it is also 
conceivable that those further away from the host centre were richer 
in gas and forming stars more actively since they formed. How- 
ever, in Fig. [8] below we show that the same trends exist in self- 
normalised galaxy properties, such as the ratio between the hot gas 
mass of a galaxy at a particular radius and at first crossing of 5 r2oo- 
This is incompatible with a scenario in which the trends presented 
here are the result of varying galaxy formation conditions: Galaxies 
are actually changing as they move closer to the host centre. 

Finally, we note that there is relatively little difference be- 
tween binning galaxies according to and Mj,ost at the snapshot 
of 'observation' (solid lines and bands) and according to the corre- 
sponding values at first crossing of 5 r2oo (M.^) and z = (M^ost) 
as shown by dashed lines. In the second case, the trends in Fig. |2] 
represent stacked evolutionary tracks of individual galaxies as they 
move towards the central host region, whereas our default binning 
method allows galaxies to change bins as a result of either the host 
mass or galaxy stellar mass varying with time. The 'fixed binning' 
method gives slightly larger fractions of gas-rich and star form- 
ing galaxies, as both the host masses and stellar masses generally 
increase with time. Measuring them earlier than the galaxy 'obser- 
vation' (M*) or later (host mass) therefore on average underpre- 
dicts the first while overestimating the second: in both cases galax- 
ies are assigned to bins with stronger environmental influence and 
therefore appear to be less strongly affected than their counterparts 
whose stellar mass and host mass were determined at the obser- 
vation snapshot. In any case, the fact that the difference is small 
implies that the trends must be predominantly caused by actual 
changes to the galaxy properties, rather than their bin designations. 



4 THE ORIGIN OF LARGE-SCALE TRENDS 

We have shown in the previous section that galaxies which are even 
moderately close to groups and clusters in the GIMIC simulation are 
systematically depleted of hot and cold gas, and are consequently 
less likely to be star-forming, compared to field galaxies of the 
same stellar mass. In the following, we aim to identify the physical 
processes responsible for these trends. There are two broad cate- 
gories of such mechanisms: direct interaction between the host and 
its galaxies at large radii on the one hand, and indirect effects such 



as pre-processing of galaxies in groups and/or 'overshooting' on 
the other. We will first address the importance of each of these in- 
direct effects, before investigating the possibility of direct galaxy- 
host interactions at large radii further below. 



4.1 Overshooting 

The first possible explanation is that a significant fraction of galax- 
ies in the outskirts of groups and clusters are not actually falling in 
for the first time. Following their pericentric passage (in most cases 
well inside r2oo)> galaxies on highly elliptical orbits may either per- 
manently escape the host on hyperbolic trajectories or at least reach 
their apocentre well beyond r2oo before falling in towards the cen- 
tral region once more. Despite their large distance from the host 
centre, these galaxies will likely already have experienced strong 
stripping due to ram pressure and tidal forces during their passage 
through the inner cluster regions and have thus lost a significant 
amount of their originally bound gas. As fewer galaxies will return 
to greater distances from the centre, this effect will naturally lead 
to a gradual increase in the fraction of star forming, gas rich galax- 
ies with increasing distance beyond the virial radius, without any 
actual environmental influence at large radii. 

The extent of this 'overshooting' in our galaxy sample can be 
judged from Fig.[3]in which we show the fraction of galaxies that 
have already ventured into the central r2oo as a function of host- 
centric distanc^l. There is no clear dependence on host mass: in all 
cases, the fraction of overshot galaxies decreases strongly between 
r2oo and 3 r2oo- Beyond this radius almost all galaxies are infalling 
for the first time. The trends in Fig. [2] are therefore likely to be 
substantially affected by overshooting at r < 3 r2oo- Beyond this 
point, however, the overwhelming majority of galaxies are infalling 
for the first time, which rules out overshooting as an explanation for 
trends in the far outskirts of groups and clusters. 

Fig.|4]highlights the difference in the radial trends for the hot 
gas fraction between overshot galaxies and those infalling for the 
first time; for simplicity and clarity, we show here trends from only 
two bins each in stellar mass and host mass. While overshot galax- 
ies are almost completely free of hot gas at all host-centric radii 
(yellow lines), the trend for those galaxies infalling for the first time 
(blue) is similar to that for the full sample shown in Fig. (2] 



4.2 Pre-processing 

A second indirect mechanism to explain environmental influence 
on galaxies at large distances from the host centre is that, in a uni- 
verse in which structures grow hierarchically, groups and clusters 
of galaxies are surrounded by other, smaller groups falling into 
them. Fig. |2] showed that strong environmental effects exist near 
the c entral regions of e v en sm a ll groups with M ~ lO'-^Mp (see 
also iBalogh & McGee' (20101), IWetzel. Tinker & Conrovl j2012h 



and|Rasmussen et al. (2014) who find a similar result in observa- 
tional data). It is possible that many of the gas-poor, passive galax- 
ies in the outskirts of a big cluster are really sitting in the central 
region of small groups, which are primarily responsible for the re- 
moval of gas. This 'pre-processing' effect has been considered as 
an explanation for reduce d star forming fract i ons beyond rono b 



a number of authors (e.g.,lBerrier et 



ILu et al.ll2012l : IWetzel. Tinker & Conrovll2012l) 



jn g tract i ons beyond rono py 
all 20091: iMcGee et aljHoogI : 



^ Note that we take into account that r2oo evolves with redshift, as the host 
grows. 
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Figure 2. Fraction of galaxies with large hot (top row) and cold gas content (middle row), and high specific star formation rates (bottom row) respectively 
at z ~ 0. Panels are split by stellar mass as given at the top; differently coloured lines represent hosts of different mass as indicated in the figure. The shaded 
regions represent statistical IfJ uncertainties. Dotted lines show the results obtained when galaxies are binned by stellar mass at infall (when crossing 5 r2oo for 
the first time), as opposed to our standard method where we bin by stellar mass at the point when the galaxy is 'observed' — there is little difference between 
the two. There are clear trends towards less gas and lower star forming Ukelihood well beyond r2oo- 



Apart from the obvious case of a galaxy actually being iden- 
tified as a satellite in an infalling group at the point of observation, 
there are two other circumstances in which a galaxy can be affected 
by pre-processing: firstly, a galaxy may have been such a satellite 
in the past, but subsequently ceased to be (for example, because it 
escaped the group's gravitational attraction or the group itself was 
tidally disrupted). On the other hand, as explained in Section 2, it 
is also possible for groups to be accreted onto a massive FoF halo 
of a galaxy cluster without being disrupted. In this case, the group 
satellite galaxies become — formally — satellites of the cluster in- 
stead of the group, which in reality continues to affect them. To 
identify such 'hidden' groups we make use of the fact that, while 
no longer forming their own FoF halo, they are nevertheless a grav- 
itationally self-bound entity which is detected as a single subhalo 
by SUBFIND, with the individual galaxies within the hidden group 
identified as sub-sub-haloes. 

As an illustration. Fig. [5] shows a projected map of all galax- 
ies within 5 r2oo from the centre of the massive galaxy cluster at 
redshift z = 0. Galaxies unaffected by pre-processing (i.e., those 
that have never been satellites in a group with M > IO'^Mq) are 
shown as black circles, filled ones representing galaxies identified 
as part of the cluster FoF halo while open ones are not. A handful 
of 'open' groups (not part of the cluster FoF halo) are shown by 
large open circles in different shades of blue, all of them at rela- 
tively large distances from the cluster centre (r/r2oo > 3). Further 
in, satellites of three hidden groups are shown by large filled cir- 



cles in shades of green. We note that, looking ahead to Fig.fTO] the 
location of these hidden groups coincides with strong overdensi- 
ties of infalling gas (as indicated by their long tails pointing away 
from the cluster centre). This confirms that these hidden groups are 
real, physical structures. Furthermore, Fig. |5] also highlights those 
galaxies which are not member of any infalling group at z = 0, but 
were in the past, with orange circles (filled/open denoting galaxies 
within/outside the cluster's FoF halo). A large number of galaxies 
in the central region (r/r2oo < 2) belong to this category, but also a 
noticeable number of galaxies around the groups at large radii. 

Combining all these various types of pre-processed galaxies, 
Fig.[6]shows that on average approximately half the galaxies around 
the massive cluster are affected, increasing towards the centre from 
~ 30 per cent at 5 r2oo to ~ 65 per cent within r2oo- In lower mass 
hosts, pre-processing is less common, decreasing to less than 20 
per cent in the case of low-mass groups. Even in this case, pre- 
processing is much more common than overshooting beyond ~ 3 
r200 and is therefore a significant contributor to the radial trends 
seen in Fig.|2] 

We finally note that we define a galaxy as having been 
been pre-processed if it was/is a subhalo of infalling group with 
Mbound > IO'-'M,!!, which in practice typically means that it has 
been within 1.5 r2oo of the group centre. Fig. |2] shows that the 
effect of environment around low-mass groups actually extends 
out slightly further than this, typically to ~ 2 r2oo- Thus, we have 
adopted a 'strong' definition of pre-processing in the present work. 
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Figure 3. Fraction of galaxies which have been within r2oo(;) of the host 
centre and therefore, if found at r > r2oo, have already passed their first 
pericentre. This is very common within 2 r2oo- Outside 3 r200: on the other 
hand, almost all galaxies are infalHng for the first time. 
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Figure 4. Radial trends in the fraction of galaxies with hot gas mass fraction 
above the threshold of 0.1, using different subsets of our infall galaxy sam- 
ple. Black lines show all galaxies, as used in Fig.|2] while blue and green 
lines show only those galaxies that are infalling for the first time and have 
also never been in a group, respectively. The green line therefore represents 
a 'clean' sample of galaxies in which any trends are due to direct interaction 
with the cluster at radii equal to or greater than the current position of the 
galaxy. Furthermore, yellow and red lines show galaxies that have already 
travelled through the inner regions of their host cluster or have at some point 
in the past been satellites in a different subhalo. Where lines do not cover 
the entire radial range this is due to a lack of galaxies in a given category 
at certain radii, e.g. overshot massive galaxies in low-mass groups at large 
radii. 
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Figure 5. Pre-processed galaxies in the massive cluster at z = (coloured 
circles). Open circles in shades of blue represent galaxies in groups identi- 
fied as separate FoF halo, while filled circles in green shades show galaxies 
in 'hidden' groups, which are part of the main cluster FoF halo, but form 
a self-bound sub-halo with total bound mass Mbound > lO'^My,. Galaxies 
shown with the same colour belong to the same group. Black circles show 
cluster galaxies which are not in any group: filled ones are identified as part 
of the cluster FoF halo, open ones reside in separate haloes. In the same 
way, orange circles show galaxies that belonged to a group only at z > 0. 
The dotted circles represent distances of (1, 2, 3, 4, 5) r2oo from the cluster 
centre. 



In principle, one could also adopt an even lower mass threshold for 
pre-processing sub-groups, but we have found that this has little 
effect on our results. 

The fraction of pre-processed galaxies which are hot gas rich 
(Mgas/M, > 0. 1) is shown in Fig.|4]together with the effect of over- 
shooting discussed above. Pre-processed galaxies are represented 
by red lines; it is clear that their hot gas fraction is significantly 
lower than in the full sample (black lines). In the case of low-mass 
galaxies (top row) the effect is virtually the same as that of over- 
shooting: in both cases, hardly any affected galaxies have signifi- 
cant amounts of hot gas, irrespective of their distance from the host 
centre. The effect is somewhat milder for massive galaxies (bottom 
row) but still quite significant. 

Combining the effects of overshooting and pre-processing, we 
can form a 'clean' sample of galaxies affected by neither: these are 
infalling for the first time directly from the field. In Fig.|4]they are 
shown in green; unsurprisingly their hot gas fractions are higher 
than in any of the four other samples, at all radii. A detailed com- 
parison between the full and clean galaxy samples, for all M» and 
Mjiost is shown in Fig.]?] The shaded bands represent the full sam- 
ple (identical to Fig. whereas the solid and dotted lines show 
the corresponding fraction of galaxies in the clean sample and their 
statistical uncertainties, respectively. 

In the case of cold gas and star formation, the trends are much 
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Figure 6. Fraction of galaxies which have been satellites of a halo with 
mass above lO'^ M ;: other than the host, as a function of galaxy position at 
time of observation. Different host mass bins are shown in different colours 
as indicated at the top of the figure. Shaded bands show the con'esponding 
statistical \a uncertainties. Galaxies in more massive hosts are more likely 
to have been satellites in a group, but with the exception of the most massive 
cluster there are no clear trends with galaxy position. 



weaker in the clean sample, and for massive galaxies with 
> IO'^'M,;, there is virtually no trend remaining - i.e., at all radii 
the fraction of galaxies with appreciable cold gas or star forma- 
tion is similar to that in the field. Where trends exist (for low-mass 
galaxies), they are strongest in the case of galaxies near the massive 
cluster. However, a much stronger environmental influence remains 
in the case of hot gas, with galaxies of all masses that we have ex- 
plored being affected out to large radii from the centres of both 
group and cluster hosts. 



4.2.1 Summary 

Our conclusions so far from this section may be summarised as 
follows: the strong radial trends in cold gas and star forming frac- 
tion seen in Fig. |2] are largely caused by overshooting and pre- 
processing, especially in the case of massive galaxies with M* 
> 10'" M0. Out of these two indirect mechanisms, overshooting is 
generally dominant within ~ 2 r200! while at larger radii the trends 
are mostly due to pre-processing. The lowest-mass galaxies, on the 
other hand, show appreciable radial trends in their retention of cold 
gas and star forming activity even when pre-processing and over- 
shooting are excluded. In terms of hot gas, the radial variation for 
both low- and high-mass galaxies is very similar in the full and 
'clean' samples. Therefore, a direct influence of the group or cluster 
environment must extend out to at least ~ 5 r2oo- In the remainder 
of the paper, we investigate this influence, and therefore use only 
the clean galaxy sample from here on. 



4.3 Direct galaxy-host interaction: tidal and ram pressure 
stripping 

In the previous section, we showed that neither pre-processing nor 
overshooting can can fully account for the depletion of cold gas and 



quenching of star formation in low-mass galaxies found in the out- 
skirts of groups and clusters. Nor can it account for the removal of 
hot gas of both low-mass and high-mass galaxies. This means that 
a process must exist by which galaxies can be influenced directly 
by the hosts at large distances from their centre. In this section, we 
argue that this process is most likely direct ram pressure stripping 
due to interaction of the galaxies with an extended hot gas halo 
surrounding groups and clusters. 

Besides ram pressure, there are several other mechanisms 
which could also be responsible for removing gas from galaxies, 
in particular tidal stripping due to the cluster potential or galaxy- 
galaxy interactions. In contrast to ram pressure, these processes 
can be expected to affect not only the (hot) gas content of galax- 
ies, but also the similarly extended dark matter haloes. In Fig.[8]we 
show the evolution of the hot gas and dark matter content of 'clean' 
galaxies falling into both the massive cluster and low-mass groups, 
focusing on low-mass (log M*/M,;, = [9.0, 9.5]) galaxies for which 
the environmental effect is strongest (see Fig.|7j, although we have 
verified that similar trends are also seen in more massive galaxies. 
Making use of our galaxy tracing results, we normalise the hot gas 
and dark matter masses by their respective values at first crossing 
of 5 r2oo- Any deviation from unity in these 'self-normalised' val- 
ues is necessarily the result of changes occuring within individual 
galaxies, and not due to potentially differing galaxy formation con- 
ditions at different distances from the host centre. 

Hot gas and dark matter clearly evolve very differently: from 
5 r2oo onwards, the hot gas mass decreases steadily with decreasing 
radius (solid lines) with an overall stronger effect in the case of the 
cluster than low-mass groups (red and black lines, respectively). 
In the former, the majority of galaxies have lost all of their hot gas 
around 3 r2oo ™d even in low-mass groups the median hot gas mass 
of galaxies is reduced to ~ 20 per cent by this point, compared to 
the value at 5 r2oo- The mass of the dark matter halo, on the other 
hand, remains nearly constant (dashed lines in Fig.[8) at these large 
radii, independent of host mass. This implies that the removal of hot 
gas in the group and cluster outskirts is due to a process targeting 
exclusively baryons while leaving the dark matter halo basically 
untouched — precisely the behaviour that would be expected from 
ram pressure stripping. 

4. 3. 1 Expected effect of ram pressure 

To test the ram pressure stripping hypothesis further, we directly 
compare the ram pressure and gravitational restoring forces on 
galactic hot and cold gas in Fig.|9] For each galaxy, ram pressure is 
computed as 

^ram = VfcM PiCM (1) 

where vjcm is the velocity of the galaxy relative to the surround- 
ing ICM of density PicM- To determine these two values, we select 
for each galaxy the N = 3000 closest gas particle^ which are not 
members of any gravitationally bound subhalo (except for the main 
subhalo in a host group or cluster). This ensures that our measure- 
ments of ICM density and velocity are not influenced by particles 
in nearby galaxies, but actually represent the ICM. To exclude con- 
tamination by gas accreted by, or stripped from, the galaxy under 

' We have experimented with various values of N and found N = 3000 to 
be the optimal value. For too low N, particle-to-particle scatter in velocity 
and density becomes noticeable, whereas at too high values of N we are no 
longer determining the local properties. 
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Figure 7. Compaiison between radial trends for the 'clean' galaxy sample, containing only galaxies which have never been within r2oo> and have never been 
satellites in another halo than the main host (shown by solid lines, dotted lines give statistical 1 a uncertainties) and the full sample from Fig.|2](shaded regions 
showing statistical Iff uncertainties, for clarity we have omitted the main trend for this sample). In the clean galaxy sample, radial trends are significantly 
weaker and start at smaller radii, with the exception of hot gas and low-mass galaxies. 



consideration, we also exclude any gas particles that have previ- 
ously been, or will subsequently be, bound to it. 

The left column of Fig.|9]shows, for 'clean' (not pre-processed 
or overshot) galaxies in low-mass groups (top) and the massive 
cluster (bottom), the distribution of resulting ram pressure values 
with varying distance from the host centre. The median trend is 
given by the thick black line, while the dark and light grey regions 
enclose 50 and 90 per cent of all galaxies, respectively. In both low- 
mass groups and the massive cluster, ram pressure is increasing to- 
wards the centre, but the trend is stronger in the latter case where 
it varies by approximately 3 orders of magnitude between 5 r2oo 
and r2oo as opposed to 'only' 2 orders of magnitude in low-mass 
groups over the same radial range. While the ram pressure expe- 
rienced by galaxies in the outskirts of both groups and clusters is 
similar (at same r/r2oo), galaxies near the centre of a cluster there- 
fore experience considerably higher ram pressure levels than their 
group counterparts. We show below that this is primarily a conse- 
quence of the higher orbital velocities of galaxies in massive clus- 
ters. Apart from this overall trend, there is also substantial scatter in 
the ram pressure values, in particular in the outer regions. Galaxies 
at a distance of 4 - 5 r2oo from the centre of a massive cluster can 
experience ram pressure differing by about 5 orders of magnitude, 
a range considerably larger than the systematic variation with ra- 
dial distance. We will investigate the origin and implications of this 
scatter in Section|5]below. 

Gas will be removed from the infalling galaxies if the 
ram pressure exceeds the gravitational restoring pressure (restor- 



ing force per unit area) exerted by the galaxy. Following 
iMcCarthv et alj bOOSi) . we compute this quantity as 

We(.)= "^^^<-^P(-^ (2) 
r 

where M(< r) is the total mass within galacto-centric radius r, p[r) 
the density of the gas phase (hot or cold, defined as discussed in 
section [3j under considerat i on, an d a is a geometric constant of 
order unitv. lMcCarthv et al] ( |2008|) find a = 2, which we adopt for 
our calculations as well, although the exact choice of this param- 
eter has no influence on our conclusions. Using equation l|2j, we 
compute hot and cold gas restoring pressure profiles for all GIMIC 
field galaxies (centrals that have never been within 5 r2oo of a group 
or cluster) as these represent the 'initial condition' of galaxies be- 
fore infall into a group or cluster. To connect the pressure compar- 
ison directly to the gas content, we furthermore find for each field 
galaxy in our sample the restoring pressure at the radius enclosing 
a series of specific hot and cold gas masses (i.e., Mgas/M») in the 
range —2.5 < logigMgas/M* < 0.5. The resulting trends, median- 
stacked in bins of similar stellar mass are shown in the middle and 
right columns of Fig.[9]and give the typical level of ram pressure re- 
quired to strip a galaxy to a given gas mass or outer limiting radius, 
respectively. We note that, in stacking galaxies, we include at each 
point only those that actually have hot or cold gas of this mass or 
extending out to this radius and only show those data points where 
this is the case for at least 5 per cent of the galaxies to give a mean- 
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Figure 8. Evolution of the hot gas and dark matter content of low-mass 
galaxies falling into hosts for the first time without having been affected by 
pre-processing. Solid lines show the median hot gas mass, dashed lines the 
median dark matter mass, both normalised for each galaxy to the respective 
values at first crossing of 5 r2oo- The red lines represent galaxies falling 
into the massive cluster in the +2a simulation, black ones those falling into 
low-mass groups. There is a clear difference between the evolution of the 
hot gas content, which decreases within ~ 5 r2oo. and dark matter, which, 
irrespective of halo mass, increases slightly until ~ 2 r2oo. This implies that 
at large distances from the host centre, the hot gas is removed by a process 
targeting exclusively baryons, such as ram-pressure stripping, and not by a 
more indiscriminate one such as tidal stripping. 

ingful picture of how tightly bound the gas typically is. For hot gas, 
we show the restoring pressure profiles over a radial range from 
to 500 kpc (physical), but because the cold gas component is much 
more centrally concentrated, we use a smaller radial range from 
to 30 kpc here. 

It is evident that cold gas is much more tightly bound than 
hot gas when comparing the middle/right top and bottom panels. 
The typical restoring pressure on cold gas ranges from 10^" Pa in 
the most massive galaxies (red) to 10^ '"^ Pa in low-mass galaxies 
(blue) at the radius enclosing 0. 1 M* in cold gas. By comparison, 
even in massive galaxies, the corresponding restoring pressure on 
hot gas is only 10^ ''^ Pa, which drops to 10^'^ Pa in low-mass 
galaxies. The reason for this is that cold gas is not only denser, but 
also sits much closer to the galactic centre. 

By comparison, the typical ram pressure reaches a maximum 
level of 10^'^'^ Pa, in the case of the massive cluster near r20Oi 
with some galaxies reaching levels up to 10^'^ Pa. This is clearly 
too low to strip cold gas in massive galaxies, but just sufficient for 
those with lower stellar masses. Outside ~ 2 r2oo and in less mas- 
sive hosts, however, no galaxies experience sufficient ram pressure 
to directly strip cold gas. Hot gas on the other hand, bound by ap- 
proximately two orders of magnitude less tightly, can be stripped 
efficiently: even massive galaxies (red) can be affected out to '--^ 2 
- 3 r2oo in clusters, and many low-mass galaxies are subject to suf- 
ficient ram pressure (~ 10^'^ Pa) even at 5 r2oo- Even in low-mass 
groups, hot gas can be expected to be stripped out to ~ 2 r2oo in all 
galaxies, and in a large fraction (> 25 per cent) out to 5 r2oo- 

These expectations agree well with the actual evolution of the 
gas content as seen in Fig. |7] In particular, stripping of hot gas 



of galaxies (top) is seen out as far as 5 r2oo in all environments, 
whereas cold gas is only affected in low-mass cluster galaxies, as 
predicted from our pressure comparison. 

We finally note that the restoring pressure profiles themselves 
show an interesting difference between cold and hot gas: those for 
cold gas are relatively flat outside ^ 5 kpc while the hot gas pro- 
files show a steady decline from the central region outwards. This 
suggests that hot gas is stripped gradually from the outside as the 
ram pressure acting on a galaxy increases, whereas when the cold 
gas finally begins to be stripped, virtually all of it will be removed 
over a short time scale. 



5 INFLUENCE OF FILAMENTS 
5.1 Origin of the ram pressure scatter 

While Fig. |9] confirms that there is a general trend to higher ram 
pressure values towards the host centre, it also reveals strong scat- 
ter, particularly in the outer regions. At r ~ 5 r2oo from the centre of 
the big cluster, the ram pressure varies between galaxies at the same 
distance from the centre by five orders of magnitude, substantially 
more than the variation in the median ram pressure over the radial 
range considered here. In this section, we investigate the origin of 
this scatter and its implications. 

An obvious possibility is that we have so far only distin- 
guished between galaxies by their radial distance from the host 
centre, thereby implicitly assuming that our hosts are spherically 
symmetric systems. This is rather unlikely: it is a long-standing 
prediction of cosmological simulations that groups and clusters of 
galaxies are triaxial systems linked by filaments of both dark mat- 
ter and gas, and there is now increasing observational evidence tha t 
this is indeed the case (e.g.. lDietrich et al.ll2012l : lAde et al.ll2012h . 
Galaxies falling in along these filaments have a very different in- 
fall experience from those accreted through largely empty regions 
(voids), which we now investigate in detail. As before, we focus 
exclusively on the 'clean' galaxy sample which are infalling for the 
first time without having been affected by pre-processing. 

As a parameter to distinguish between galaxies in filaments 
and voids, we choose the 'local overdensity', which we define as 

= Plocal / Pprofile (''galaxy ) (3) 

where piocal is the locally determined ICM density obtained as dis- 
cussed in the previous section and Pprofile (''galaxy) is the correspond- 
ing spherically averaged gas density at the radius of the galaxy. In 
a perfectly spherically symmetric host, each galaxy would have a 
value of Ap = 1 ; in a more realistic system, galaxies in filaments 
are those with the highest Ap while those in voids have the lowest 
values. Having determined Ap, we then bin galaxies in each snap- 
shot according to host mass and rank the galaxies in each bin by Ap . 
The highest quartile in each bin is identified as filament galaxies, 
and the lowest as those in voids. This ensures that we have equal 
numbers of filament and void galaxy data points and that both also 
have the same distribution in redshift and host mass. In Fig.llOlwe 
show the locations of both filament and void galaxies in the mas- 
sive cluster at z = 0, superimposed on a map of the ICM gas density. 
As expected, filament galaxies are strongly spatially clustered and 
are mostly found in two bands fanning out towards the top right 
and the left, which is exactly the region where prominent filaments 
can be seen in the gas density map. The void galaxies, however, are 
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Figure 9. Comparison of ram and restoring pressure on galactic gas. For ram pressure (left), the top panel shows low-mass groups, the bottom panel a massive 
cluster. The solid black line shows the median trend, dai'k and light grey bands enclose 50 and 90 per cent of galaxies. For restoiing pressure (middle/right 
columns), the top panels show hot gas whereas the bottom panels show cold gas; in both cases differently coloured hues show galaxies of different stellar 
mass. In the middle column, the restoring pressure is plotted as a radial profile, in the right as a function of enclosed specific gas mass. See text for details. 
This comparison predicts substantial stripping of hot gas in hosts of all masses. Cold gas, on the other hand, is only expected to be significantly stripped in the 
central cluster regions. This agrees well with what is observed in the GIMIC galaxies. 



much more evenly spread and tend to avoid the filament regions^ 
We take this as confirmation that our classification based on Ap 
is a physically meaningful distinction between void and filament 
galaxies. 

As a first step, we show in Fig.[TT]a decomposition of the ram 
pressure values from Fig. [5] (right column) into ICM density (left) 
and its velocity relative to the galaxy (middle). Galaxies near low- 
mass groups are shown on top and those around the massive cluster 
in the bottom row. In both cases, the median trends for filament 
galaxies are represented by red lines, void galaxies in blue, while 
the green solid line shows the overall trend (irrespective of Ap). 
In the left column showing ICM density, we furthermore show the 
median-stacked spherical density profile as a dotted green line; in 
the middle column the velocity of the galaxy relative to the host 
centre-of-mass frame is shown in black. 

The difference in ICM density between filament and void 
galaxies is evidently quite strong with the latter around an order 
of magnitude lower. The difference is least pronounced near r2oo 
(0.6 and 0.9 dex in groups and the cluster, respectively), and in- 
creases outwards, reaching 1.1 and 1.5 dex at 5 r2oo- In both cases, 
the density trends are relatively flat, in particular in the case of fil- 
ament cluster galaxies where there is no significant radial variation 

Note that this is a projected map, so that some blue points (void galaxies) 
may still appear to lie in regions of high gas density. 



outside r2oo- The increased density difference between voids and 
filaments towards larger radii is driven mostly by a moderate de- 
crease in the density around void galaxies. Including scatter within 
each band, ICM densities around 5 r2oo vary by approximately 1.5 
orders of magnitude in groups and two in clusters, thus account- 
ing for about half of the overall ram pressure scatter. As expected, 
the typical ICM density for all galaxies (green) lies approximately 
halfway between the filament and void trends. It is worth pointing 
out that the overall spherical density profile (green dotted) is gen- 
erally lower than the ICM near galaxies, indicating that galaxies 
preferentially live in overdense regions. This is particularly true in 
low-mass groups, where the spherical density is even lower than 
the median density around void galaxies — in this environment, 
'voids' really represent the average density, whereas filaments are 
highly overdense regions. In clusters, on the other hand, voids are 
actually underdense with respect to the spherical average. 

The ICM velocity behaves in a rather different way. While the 
overall radial variation in low-mass groups (middle top panel, green 
line) is comparable to that in density, it is much stronger in the case 
of the cluster (middle bottom panel), where it varies by significantly 
more than an order or magnitude, decreasing with increasing ra- 
dial distance. The velocity between a galaxy and its surrounding 
ICM is almost always considerably lower than that with respect to 
the host centre-of-mass frame (black), the difference being great- 
est at largest radii (more than a factor of 10 at its most extreme, ~ 
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Figure 10. Gas filaments around the massive cluster at z = (background 
image, lighter colours corresponding to higher gas surface density). There 
is a good correlation between dense regions and the location of galaxies that 
we identify as in filaments (red points). Void galaxies (blue points) show no 
such strong clusteiing. The dotted white lings mark projected radii of 1 - 5 
1200- 



1000 and 100 km/s, respectively). The closest match between both 
is found for void galaxies near r2oo where the two velocities virtu- 
ally agree, indicating that in this case the ICM is largely stationary 
in the host frame, whereas it is infalling together with the galaxies 
at larger radii. Perhaps surprisingly, the radial trends are stronger 
for void galaxies, so much so that at 5 r2oo. their velocity relative to 
the ICM is actually lower than in filament galaxies at the same ra- 
dius outside ~ 3 r2oo- But even within this radius, the much higher 
ICM density around filament galaxies means that they experience 
higher levels of ram pressure (see right column) throughout their 
infall than void galaxies, with no evidence of 'dynamical shield- 
ing' by which the co-flow of filament gas could protect galaxies 
from ram pressure by the surrounding gas (except possibly very 
close to r2oo> see below). 

As the velocity enters into ram pressure in quadrature, its 
strong radial trend is the main cause for the radial variation in ram 
pressure (accounting for a variation of ~ 1 order of magnitude 
in Pram for filament galaxies and three orders for those in voids). 
Looking at the right-hand column, almost the entire extent of radial 
variation in ram pressure can therefore be attributed to an increased 
velocity nearer the centre, while its scatter is strongly influenced 
by variation in ICM density at the same radius. We next explore the 
implications of this on the properties of void and filament galaxies. 



5.2 Effect on gas fractions 

The results above show that galaxies in filaments will, on average, 
experience higher levels of ram pressure than those in voids. How- 
ever, this does not necessarily imply that these galaxies are also less 
likely to contain significant amounts of gas: it is plausible, for ex- 
ample, that galaxies falling in along a filament could accrete more 
gas prior to stripping, since far away from the host centre gas is 



evidently much more abundant within filaments than outside (see 
Fig.[il. 

To see which influence is stronger, we plot in Fig.[T2]the frac- 
tion of galaxies with substantial hot gas (logio Mj,ot (.a.s/M» > —0.5) 
as a function of distance from the host centre^ Again, we distin- 
guish between galaxies in filaments and voids (red and blue respec- 
tively, with green lines showing the overall sample), for both low 
mass groups (top) and the massive cluster (bottom). The solid lines 
show the fractions themselves, while the shaded bands represent 
the statistical (Poisson) uncertainties. As we are now once more 
looking at internal galaxy properties, we split our sample by stellar 
mass, each bin represented in a different column. 

The clearest trends are visible for low mass galaxies with logio 
(M^/Mq) < 10.0. In both the big cluster and low-mass groups, they 
reveal a clear influence of local environment, in particular in the 
outer regions (beyond 3 r2oo)> in the sense that galaxies in fil- 
aments (red) are much less likely to be hot gas rich than those in 
voids (blue): even at a distance of 5 r2oo from the cluster centre, 
only ~ 5 per cent of filament galaxies in the lowest stellar mass 
bin are hot gas rich, compared to ~ 50 per cent of void galaxies 
(bottom left panel). Evidently, the increased level of ram pressure 
stripping experienced in filaments is more than enough to negate 
any potential benefit in gas accretion these galaxies might have ex- 
perienced. The difference is slightly smaller in low-mass groups, 
but even here a high hot gas content is found in only ~ 10 per cent 
of filament galaxies as opposed to ~ 50 per cent in voids. In the 
case of low-mass galaxies, the influence of filaments is therefore 
clearly an increased stripping of hot gas at large radii. 

Note that like all results presented in this paper, these trends 
are found at low redshift (z < 0.5). We have verified that at high 
redshift {z = 2), an equally large fraction of our filament and void 
galaxies are hot gas rich. The clear discrepancy even at 5 r2oo be- 
tween the two galaxy populations at low redshift therefore indi- 
cates that the stripping influence of filaments extends even beyond 
5 r2oo> and also implies that this effect is predominantly a low- 
redshift phenomenon. We intend to pursue this result in a future 
study, together with the general redshift evolution of environmen- 
tal trends which is now increasingly probed observationa lly (e.g., 
iBalogh et alj201ll : lMcGee et al.ll201ll : ISnvder et al.ll2012h . 

For galaxies of higher stellar mass, the smaller sample size 
makes it more difficult to draw statistically robust conclusions as 
to the difference between filament and void galaxies. Any potential 
difference is certainly smaller than in the lower stellar mass bins — 
within the large statistical uncertainties, almost all trends are con- 
sistent with no difference at all between filament and void galaxies. 
It is particularly unfortunate in this respect that there are very few 
high-mass void galaxies: in the same way as galaxies in general are 
preferentially found in overdense regions (see above), the filaments 
are clearly a more attractive habitat for massive galaxies than voids, 
as is evident from the different magnitude of the errors. There may 
be a hint of evidence for an increased hot gas fraction in filament 
galaxies near r2oo, in particular in the massive cluster (bottom row). 
Keeping in mind that the ram pressure levels between filament and 
void galaxies are most closely matched at small radii, this might 
therefore be an indication for a small amount of dynamical shield- 
ing by filaments once galaxies approach r2oo- In any case, this dif- 



^ This threshold is a factor of ~ 3 higher than that which we have used 
above. While both limits give rise to similar overall trends, we have found 
that this threshold value leads to the clearest separation between void and 
filament galaxies. 
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Figure 11. Density (left column), relative velocity {middle column) and resulting ram pressure (right column) of the ICM around 'clean' galaxies. Low-mass 
groups are shown at the top, massive clusters in the bottom row, each split into filament galaxies (red) and those in voids (blue) as described in the text; green 
lines represent the undivided sample. Also shown is the stacked spherical density profile (green dotted) and the velocity with respect to the host centre of mass 
frame (black). The median and 50 % spread of each bin are shown by solid lines and shaded regions respectively; the darker shaded region around the median 
line gives the statistical uncertainty in the median (which is small due to the large sample size). The radial trends in ram pressure are mostly due to radial 
variation of velocity, whereas the scatter is dominated by variation in local density. 



ference is small compared to the very strong influence of filaments 
in the outer regions. 



6 SUMMARY AND DISCUSSION 

We have analysed the gas content and star formation rate of galax- 
ies in the vicinity of groups and clusters in the GIMIC suite of 
hydrodynamic cosmological simulations. By using zoomed initial 
conditions, these simulations include rare objects such as a massive 
galaxy cluster while simulateneously resolving individual galaxies. 
We have focussed on determining how far out from the host centre 
galaxies are environmentally affected, as well as finding the phys- 
ical mechanisms behind this influence. Our main conclusions may 
be summarised as follows: 

(i) The fraction of galaxies with large hot and cold gas fractions 
(Mgas/M» > 0.1) decreases systematically with decreasing distance 
from the host centre; the same is true for the fraction of star form- 
ing galaxies (sSFR > 10^^' yr^^). These trends extend out to very 
large distances from the host centre (up to ~ 5 r200i corresponding 
to ~ 10 Mpc for a massive cluster). 

(ii) In terms of cold gas and star formation, the radial trends 
are explained by galaxies past first pericentre ('overshooting', es- 
pecially within ~ 2 — 3 r2()o) and those having been 'pre-processed' 
in infalling groups (~ 10 — 50 per cent of galaxies, increasing with 
increasing halo mass. However, even amongst not pre-processed 
galaxies infalling for the first time radial trends extending out to 5 



r2oo are still seen for the lowest-mass galaxies we have examined 
(9.O<logi,)(M*/M0) <9.5). 

(iii) In the case of hot gas, the radial trends cannot be explained 
solely by overshooting and/or pre-processing. This implies a direct 
interaction with the host group/cluster out to very large radii. 

(iv) We have shown that the ram pressure exerted by the ICM 
is strong enough to significantly strip the hot gas haloes around 
both low- and high-mass galaxies out to several multiples of r2oo- 
However, it is not strong enough to significantly strip cold gas far 
beyond r2oo (except for low-mass galaxies with M» < 10^'^ Mq in 
a massive cluster). This agrees well with the actual gas content of 
GIMIC galaxies. 

(v) At large distances from the host centre (r > 2 r2oo) there is 
a substantial difference between the velocity of galaxies relative 
to the host centre-of-mass frame and relative to their surrounding 
ICM (the latter being lower), increasing with host mass and dis- 
tance from the host centre. For low-mass groups (M^ost ~ lO'-^ Mq) 
the difference is a factor of ^ 3, increasing to ~ 10 for the massive 
cluster (Mjiost ~ lO'^ Mq) simulated in GIMIC. This implies that 
calculations based on the centre-of-mass velocity will significantly 
overestimate the ram pressure exerted on galaxies at large r/r20()- 

(vi) At large distances from the host centre, there is very large 
scatter in the ram pressure values experienced by different galaxies 
at similar rk2oo, reaching a factor of ~ 10^ in the case of the mas- 
sive cluster. This scatter is dominated by variations in ICM density: 
in filaments, ram pressure can be up to two orders of magnitude 
larger than in the lowest-density regions. 

(vii) Because they experience higher ram pressure, low-mass 
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Figure 12. Radial trends in the fraction of liot gas rich (logio Mhotsas/M* > —0.5) 'clean' galaxies of varying stellar mass (increasing left to light) in small 
groups (top) and the massive cluster (bottom). Red lines represent filament galaxies, blue those in voids, and green all galaxies. For low-mass galaxies in 
particular, galaxies in filaments are significantly less likely to have a high hot gas content than those falling in through voids. For more massive galaxies, 
however, there is no significant difference between the two infall routes. 



galaxies (logio M,/M0<9.5) in filaments are significantly 
less likely to contain substantial amounts of hot gas (logio 
Mhot gas/M* > —0.5) than those infalling through voids at the same 
distance from the host centre (~ 5 vs. ~ 50 per cent at a distance 
of 4 - 5 r2oo from the massive cluster for the lowest mass galaxies 
in our sample). The difference is less pronounced for more massive 
galaxies. 

Our result that the properties of galaxies near groups and 
clusters but outside r2oo differ from those at much greater dis- 
tances agrees with observational studies, which are increasingly 
supporting a picture in which t he influence of envi r onment reaches 
well beyond this radiu s (e.g., lBaloghetal.lll999l ; ILu et al.|[201^ ; 
iRasmussen et alj|2012h . So far, pre-processing and overshooting 
have been advocated as th e most likely explanations for thes e trends 
(e.g., lBaloghetalJll999l : IWetzel. Tinker & ConroyI |2012|) which 
agrees with our GIMIC results in the case of cold gas and star for- 
mation activity. However, the simulations also predict that a direct 
interaction with the host and its filaments is sufficiently strong to 
remove the hot gas of galaxies as far out as 5 r2oo- Testing this ob- 
servationally in a direct fashion will be challenging, though, as it is 
extremely difficult to probe the bulk of the hot baryons around even 
very local galaxies. 

We finally note that, even though we conclude that it is un- 
likely for cold gas to be stripped from galaxies outside r2oo> this 
does not mean that it remains wholly unaffected: the removal of hot 
gas stops the replenishment of cold gas lost through star formation, 
and will therefore have some indirect impact on the cold gas con- 
tent and star f ormation activity of gal axies, although not immedi- 
ately (see also lBlanton et alj200d and lWilman, Zibetti & Budavaril 
l201Ct who show that the fraction of red galaxies depends largely 
on the local, not the large-scale galaxy density). Likewise, a de- 
crease in star formation will, after a further delay, leave an imprint 



in the galaxy colour. As well as accurately modelling the removal 
of hot gas inside r2oo. taking into consideration the loss of gas be- 
fore even reaching this radius may further improve the accuracy of 
semi-analytic models of galaxy evolution. 

As galaxies move inside r2oo and fall deep within the poten- 
tial wells of their hosts, they experience not only increased levels 
of ram pressure, but also other effects such as tidal stripping and 
harassment by other galaxies. We have seen in the case of overshot 
galaxies that these effects have a dramatic impact on their internal 
properties. In a future paper we explore in detail this complex in- 
terplay of environmental influences in the central region of groups 
and clusters. 
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